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■ ABSTRACT 
Oh! 

. Context. Spectral characterization of sub-stellar companions is essential to understand their composition and formation processes. 

■ However, the large contrast ratio of the brightness of each object to that of its parent star limits our ability to extract a clean spectrum, 
free from any significant contribution from the star. During the development of the long slit spectroscopy (LSS) mode of IRDIS, the 
dual-band imager and spectrograph of SPHERE for the Very Large Telescope (VLT), we proposed a data analysis method to estimate 

. and remove the contributions of the stellar halo and speckles to the spectra. 

Ph ■ Aims. This method has never been tested on real data because of the lack of instrumentation capable of combining adaptive optics 

PQ ' (AO), coronagraphy, and LSS. Nonetheless, a similar attenuation of the star can be obtained using a particular observing configuration 

where the slit is positioned on the faint companion while keeping the bright primary outside. 

Methods. Test data were acquired with this slit configuration using the AO-assisted spectrograph VLT/NACO. We obtained new 
/- and //-band spectra of SCR Jl 845-6357 B, a T6 companion to a nearby (3.85 + 0.02 pc) M8 star. This system is a well-suited 
Q ' benchmark as it is relatively wide (~1.0") with a modest contrast ratio (~4 mag), and a previously published JHK spectrum is 

^ \ available for reference. 

^ . Results. We demonstrate that (1) our method is efficient at estimating and removing the stellar contribution, (2) it allows to properly 

^ ■ recover the overall spectral shape of the companion spectrum, and (3) it is essential to obtain an unbiased estimation of physical 

parameters. We also show that the slit configuration associated with this method allows us to use long exposure times with high 
throughput producing high signal-to-noise ratio data. However, the signal of the companion gets over- subtracted, particularly in our 
/-band data, compelling us to use a fake companion spectrum to estimate and compensate for the loss of flux. Our spectral analysis 
provides an estimation of Te/f = 1000 + 100 K leading to R = 0.7 + 0.1 Rjup, a value that closely agrees with evolutionary models 
for ages older than 1.5 Gyr. Finally, we report a new astrometric measurement of the position of the companion (sep = 0.817", PA = 
' 227.92 deg), which has undergone a significant proper motion since the previous measurement. 

(N : 

_ Key words, stars: binaries: close - stars: brown dwarfs - astrometry - methods: data analysis - instrumentation: adaptive optics 

o 

■ 1 ■ Introduction (Macintosh et al. 2008) for Gemini, where the images and spec- 
tra of exoplanets will be strongly limited by the speckle noise 

Near-mfrared emission spectra of sub-stellar companions di- associated with the star (Soummer et al. 2007). 
rectly imaged around pre-main sequence and main sequence 

stars is crucial when it comes to characterizing their atmospheric In the context of the development of the differential spectro- 
(e.g. composition and cloud properties) and physical (effective imager SPHERE/IRDIS (Dohlen et al. 2008), we proposed in 
temperature Te/f, surface gravity log g, radii, mass and lumi- Vigan et al. (2008) a new data analysis method to estimate and 
nosity) properties, as well as understanding their formation and subtract the stellar halo and speckles in high-contrast long slit 
evolutionary processes. Obtaining high-quality spectra for faint spectroscopy (LSS) data. The method was designed to attenu- 
companions at close angular separations from bright stars is ate the speckle noise and extract a clean companion spectrum, 
challenging, as one needs to overcome the large contrast ratio free from any significant stellar contribution. It has proven to 
in the brightnesses of the two objects. Low-mass companions be very efficient in realistic simulated data for SPHERE/IRDIS, 
are usually only slightly brighter than the point spread function but has never been tested in real data. To demonstrate the 
(PSE) wings, making it necessary to develop methods to estimate method, we proposed using the VLT/NACO AO-fed spectro- 
and subtract at all wavelengths the stellar contribution, which graph (Lenzen et al. 2003; Rousset et al. 2003) in a particular 
can bias the spectral analysis. This will be even truer for up- configuration where the slit is positioned on a faint sub-stellar 
coming high-contrast coronagraphic imagers such as SPHERE companion and the bright primary is kept outside of the slit 
(Beuzit et al. 2008) for the Very Large Telescope (VLT) or GPI to avoid saturation. Although not strictly identical to what will 
be available in SPHERE/IRDIS, this configuration improves the 



* Based on observations collected at the European Southern dynamical range compared to the more standard configurations 
Observatory, Chile, ESO program 486.L-0077. where the bright primary and the sub-stellar companion are kept 
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Table 1. SCR 1845 B astrometry 



Date 


Separation 


PA 


Ref. 




(as) 


(deg) 




1 May 2005 


1.176 + 0.001 


170.22 + 0.08 


(1) 


28 May 2005 


1.170 + 0.003 


170.20 + 0.13 


(2) 


2 May 2006 


1.064 + 0.004 


177.20 + 0.06 


(3) 


20 Oct. 2010 


0.817 + 0.003 


227.92 + 0.46 





References: (1) Montagnier et al. (2006); (2) Biller et al. (2006); (3) 
Kasper et al. (2007). 



within the slit (e.g. Mohanty et al. 2007; Lafreniere et al. 2008; 
Janson et al. 2010). 

SCR J1845-6357 - hereafter SCR 1845 - is a close (3.85 ± 
0.02 pc; Henry et al. 2006) M8.5 star (Henry et al. 2004) discov- 
ered in the context of the SuperCOSMOS survey (Hambly et al. 
2001). Biller etal. (2006) discovered a bound companion or- 
biting at -4.5 AU from the primary using VLT/NACO AO- 
imaging, making it the first system with a tight brown-dwarf 
companion to a main-sequence star. They estimated SCR 1845 B 
to be a T4.5-T6.5 brown dwarf based on its colors in narrow- 
band filters centered inside and outside of the 1.6 yum CH4 ab- 
sorption band ("spectral diff'erential imaging" mode). The T6 
spectral type was later confirmed using NACO spectroscopic 
observations (Kasper et al. 2007), which allowed them to con- 
strain Teff and log g to be -950 K and 5.1 dex respectively, 
assuming solar metallicity. They also estimated an age ranging 
from 1.8 Gyr to 3.1 Gyr for this system, which translates into 
a mass of 41.0-51.8 Mjup for the companion using the evolu- 
tionary models of Burrows et al. (1997). SCR 1845 B is then 
very similar in terms of distance, age, and spectral type to the 
well-known T1/T6 dwarf binary orbiting the K4.5V star e Indi A 
(McCaughrean et al. 2004; King et al. 2010). 

Systems such as SCR 1845 and e Indi Ba/Bb are extremely 
appealing for the calibration of evolutionary models: their prox- 
imity to the Sun and their small orbital radius make them ex- 
cellent candidates for an independent, dynamical estimation of 
the system mass on relatively short timescales. It is all the 
more important that evolutionary tracks are not well-calibrated 
at low masses and young ages (e.g. Konopacky et al. 2010). In 
the present work, we provide a new astrometric measurement 
of the position of SCR 1845 B, which has a large proper mo- 
tion with respect to previous measurements, and we obtain new 
VLT/NACO /- and //-band spectra to test our data analysis 
method. We start by giving a brief description of the observa- 
tions and basic data reduction steps in Sect. 2, before describing 
in detail the results, performance, and problems of our data anal- 
ysis in Sect. 3. In Sect. 4, we discuss our new spectral extraction 
of SCR 1845 B, and finally in Sect. 5 we discuss the impact of 
suppressing the stellar contribution and compare our data anal- 
ysis method to previously published methods in the context of 
future high-contrast spectrographs. 



2. Observations and data reduction 

2.1. Imaging 

Data were obtained on 20 October 2010 with the instrument 
VLT/NACO in service mode (ESO program 486.L-0077). The 
infrared wavefront sensor with 14 x 14 lenslets array was used 
together with a dichroic transmitting 90% of the light in K band 
to the NAOS wavefront sensor and 90% of the light in either J or 
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Fig. 1. VLT/NACO image of SCR 1845 AB obtained in 7^-band 
with the S27 camera on 20 October 2010. North is up and east 
is left. The center of the primary is symbolized by a white 
star symbol. The relative positions of SCR 1845 B reported by 
Biller et al. (2006) and Kasper et al. (2007) are overlaid on the 
image, respectively, as a green plus sign and blue cross. The pe- 
culiar orientation adopted for the slit in our work is also overlaid 
on the image. 



H band to CONICA. We took advantage of the median-quality 
atmospheric conditions in the /- (seeing of -l.T', airmass of 
-1.5) and H- (seeing of -0.8'', airmass of -1.4) bands, provid- 
ing a good AO correction. 

A series of 5x5 sec exposures were recorded with the S27 
camera and the //-band near-infrared filter, providing a spatial 
scale of 27.15 mas/pixel. These frames were divided by a lamp 
flat-field exposure and corrected for the efl'ects of bad pixels 
using customized IDL routines. The cleaned frames were then 
median-combined to produce a final image (see Fig. 1). Since 
we did not foresee performing precise photometry, no dithering 
procedure was applied during the observations and no sky frame 
was subtracted. Precise relative astrometry was performed by 
fitting two-dimensional (2D) Gaussian functions to both com- 
ponents of the system. Table 1 compares our relative astrometry 
to previously reported values. Given the proximity of SCR 1845 
to the Sun and the short angular separation of the system, the 
orbital motion has been significant in the past five years. 



2.2. Spectroscopy 

NACO was used in LSS mode to obtain R = 400-500 
spectra in /- and //-band with the 86 mas slit. The spec- 
troscopic modes used were S54-4-SJ for /-band, providing 
~2 nm/pixel spectral dispersion over the range 0.91-1.40 yum 
and 54.60 mas/pixel spatial resolution, and S27-4-SH for H- 
band, providing ~ 1 nm/pixel spectral dispersion over the range 
1.37-1.84 Jim and 27.15 mas/pixel spatial resolution. The slit 
was aligned in a very particular way to ensure that the bright 
core of the PSF of the primary did not enter the slit (see Fig. 1), 
allowing a remarkable increase in sensitivity to the companion 
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Fig. 2. Comparison of the /- (left) and ^-band (right) spectra before and after the spectral deconvolution data analysis used to 
remove the contribution of the star. The top panels display the spectra with no star subtraction, the middle panels the spectra after 
data analysis at the same display scale as the top panels, and the bottom panels the spectra after data analysis with an appropriate 
display scale to identify the areas of over- subtraction (darker areas). The areas of strong over-subtaction are identified using small 
white arrows in the bottom panel. The display scales are different in the /- and //-bands. 



signal, while keeping signal from the star, which will be esti- 
mated and subtracted during the data analysis. A standard nod- 
ding sequence of 10'' was used to obtain eight exposures of 
70 sec (DITxNDIT = 10x7), resulting in a 560 sec total ex- 
posure time, in both the /- and //-bands. A telluric calibrator 
(HIP 89729, spectral type 05.5III) was observed just after the 
object at a similar airmass. A spectrum of the NACO internal 
argon arc lamp was recorded for wavelength calibration. 

The spectroscopic data was reduced using RED SPEC \ the 
data reduction software originally designed for the high- 
resolution spectrograph Keck/NIRSPEC. The software was used 
to perform the following data reduction steps: division by flat- 
field, correction of bad pixels, correction of the spectrum distor- 
tion, and wavelength calibration. The diff'erence between pairs 
of corrected frames was obtained to eliminate the sky contribu- 
tion. The frames were recentered on the same spatial position 
using the signal of the companion as a reference, and finally 
median-combined to create the final reduced spectrum. The re- 
sulting spectra in /- and //-band are visible in the top panels of 
Fig. 2. 



^ http : //www2 . keck . hawaii . edu/inst/nirspec/redspec . html 



3. Spectral deconvolution 

3.1. Context 

In high-contrast imaging and spectroscopy, the signal of the 
companion is usually close to or below the brightness level of 
the stellar halo and speckles, requiring the use of specific data 
analysis methods to alleviate the contribution of the star, which 
otherwise contaminates the spectral extraction of the companion. 
This is extremely important (i) to accurately retrieve the shape of 
the companion continuum and (ii) to isolate the spectral features 
of the companion. In Fig. 2 (top panels), we see that the com- 
panion signal (thin straight line) is overlaid on the stellar halo 
and speckles. The signal of the companion is then biased by the 
star signal, so the goal of the data analysis step will be to ob- 
tain a spectrum of the companion that is free of any significant 
contribution from the star. 

The diff'erentiation of point source spectra from a constant 
or varying background is a problem encountered in very dif- 
ferent contexts with LSS, and several solutions have been pro- 
posed. An iterative approach was proposed by Lucy & Walsh 
(2003) and successfully used to extract clean supernova spec- 
tra (Blondin et al. 2005). This method is based on a two-channel 
algorithm that restores a PSF-like component in a 2D image 
and an underlying extended background separately. However, 
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Fig. 3. Illustration of the various data analysis methods that have been proposed and used to estimate and subtract the stellar halo 
and speckles at the position of a faint companion. The star, halo, and speckles are represented in dark grey, respectively, with a 
five-pointed star, a wide radial gradient, and small grey circles; the companion is represented by a smaller red five-pointed star; the 
portion of starlight used to calibrate the halo/speckles at the position of the companion is represented by a hashed green area. The 
chromatic motion of the halo and speckle pattern is represented by dark grey arrows starting from the speckles and pointing outwards 
(i.e. motion for increasing wavelength). A short explanation of the methods is given beneath each illustration. In the Vigan et al. 
(2008) configuration, the slit is coupled with a Lyot coronagraph, represented by a circular focal-plane mask at the center of the slit. 



it requires precise knowledge of the slit- spread function (SSF) 
to restore the point-source channel, and a misevaluation of 
the width of the spatial resolution kernel necessary for the 
point- source/extended background discrimination will aff'ect the 
restoration (see Blondin et al. 2005, Sect. 3). This method has so 
far never been used to extract faint sub- stellar companion spec- 
tra. 

More specifically to sub- stellar companions. Close et al. 
(2005) used the symmetry of the stellar profile to extract the 
spectrum of AB Dor C: they subtracted the contribution of 
AB Dor A at the position of the companion using a spectrum ac- 
quired with the telescope derotator rotated by 180°. For the spec- 
trum of 2M 1207 b, Mohanty et al. (2007) modeled the stellar 
halo in each spectral channel using a sum of three Gaussians plus 
a second order polynomial. For the spectrum of IRXS 1609 b, 
Lafreniere et al. (2008) simply interpolated the star contribu- 
tion with a straight line underneath the companion signal. This 
assumption, although quite simplistic, is well-justified by the 
wide separation of the system (2.22'0. The same approach was 
used by Kasper et al. (2007) for the spectrum of SCR 1845 B. 
Figure 3 presents an illustration of these diff'erent data analysis 
methods with the relevant information: position of the star, halo 
and speckles, position of the companion, orientation of the slit, 
chromatic motion of the speckles, and portion of starlight used 
to calibrate the halo and speckles. 

At higher contrasts, we developed a method for the 
SPHERE/IRDIS LSS mode (Vigan et al. 2008) based on the 
"spectral deconvolution" - hereafter SD - originally proposed 
by Sparks & Ford (2002). It is based on the linear dependence 
of the position and size of all diff'raction-related patterns (stellar 
halo, speckles. Airy rings) on wavelength, when the position of 
an object in the field, such as a sub-stellar companion, is fixed. 



The main diff'erence from the present work is that IRDIS will 
combine LSS with coronagraphy to be able to characterize faint 
companions around nearby bright stars. In the case of IRDIS, 
the diffraction patterns will then move away radially from the 
star inside the slit with increasing wavelength. 



In the present work, we slightly tilted the slit to decenter the 
star (see Fig. 3), allowing us to use long exposure times with- 
out saturating the star signal. In this configuration, the speckles 
still move away radially from the star, but because of the tilted 
slit they do not necessarily remain inside the slit over the full 
wavelength range. This purely gometrical eff'ect is illustrated in 
Fig. 3 (rightmost illustration): the closer a speckle is to the star, 
the more orthogonal to the slit its chromatic motion is. The re- 
sult is that speckles at small angular separations will either en- 
ter or exit the slit, creating a variation in luminosity in the final 
spectrum that cannot be calibrated, while speckles at larger sep- 
arations (e.g. close to the companion) will remain in the slit over 
a wider range of wavelengths, allowing an a posteriori calibra- 
tion of their luminosity variation with wavelength. Another way 
of stating the problem is to say that because of the perspective 
eff'ect, the slit is much wider at larger angular separations from 
the direction of the star. In our data, the separation of the com- 
panion is large enough (0.817'') to ensure that this eff'ect does 
not prevent us from applying the SD. Moreover, at this separa- 
tion the data is mostly dominated by the star halo, which shows 
smoother variations than bright speckles very close to the star. In 
the following sections, we focus on validating our data analysis 
method assuming that this eff'ect does not seriously hamper the 
analysis at the position of the companion. 
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Spectral Deconvolution - ^ = 1 .65 lam Influence of over-subtraction - H band 
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Fig. 4. Comparison of cuts along the spatial dimension of the 
spectra before (plain, black) and after (plain, red) SD at /I = 
1.65 yum. Since the star was kept outside the slit to prevent sat- 
uration, the profile of the primary before SD does not show the 
true peak of the primary PSF, but instead the level of the halo 
brightness with the slit decentered by ~0.2''. The PSF profile 
from the star obtained from Fig. 1 is overplotted (dashed, black) 
to give an estimation of the true contrast (4.17 mag in //-band) 
between the companion and the primary. For clarity, the main 
plot has been cut at 16 ADU, but the sub-plot shows the full 
profiles on a log scale with the same units in x and y. 



3.2. Application to the data 

The SD was applied separately to the /- and //-band spectra 
obtained at the end of the reduction steps detailed in Sect. 2. 
The result is displayed in the middle panels of Fig. 2. Compared 
to the spectra before SD, the contribution of the star has been 
significantly attenuated, although some residuals are still visible 
at the position of the star. This is particularly visible when using 
an appropriate display scale (Fig. 2 bottom panels). Two eff'ects 
are clearly visible: 

- The non-zero residuals at the position of the star signal 
because, as previously explained, the diff'raction-related 
patterns move through the slit. This eff'ect makes it difficult 
to precisely evaluate their contribution close to the star, 
where their chromatic motion is almost orthogonal to the 
slit and they are brightest. This effect is significant within 
±0.5-0.6'' of the position of the star. 

- the over- subtraction visible on either side of the companion 
spectrum because the estimation of the stellar contribution 
at the position of the companion in the rescaled spectrum is 
biased by the companion signal. To mitigate this eff'ect (but 
not cancel it), the companion signal is masked in the fitting 
process with a 5A/D aperture. 

The residuals at the position of the star are not problematic be- 
cause we attempt to extract the spectrum of the faint companion, 
so we do not study its effect in any more detail. While in H band 
the over- subtraction appears as thin areas either below or above 
the companion spectrum (see white arrows in Fig. 2), it appears 
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Fig. 5. Comparison of the SD-only and SD-h compensation spec- 
tra in //-band. The SD-Fcompensation spectrum corresponds to 
the SD-only spectrum corrected for the over- subtraction by in- 
troducing a fake companion into the data. 



much stronger in /-band, especially at the center and the red- 
end of the spectrum, making it impossible to recover a spectrum 
over the full band (see Sect. 3.4). The over- subtraction is ana- 
lyzed in more detail in the next section. Figure 4 illustrates the 
effect of the data analysis with a cut of the spectrum before and 
after SD in //-band Sit A = 1.65 fim. Since the star was kept out- 
side the slit, the profile before SD actually shows the level of the 
star halo when the slit was decentered by ~0.2''. The level of the 
PSF profile obtained from Fig. 1 is overplotted to show the true 
brightness contrast difference of the companion with respect to 
the primary (4.17 mag in //-band). The effect of the SD is clearly 
visible, as the star contribution was significantly reduced, leav- 
ing mostly residuals with an average close to zero. On either side 
of the companion signal, the over- subtraction appears as nega- 
tive areas of 0.r'-0.2'' width. Apart from these areas, the signal 
of the companion is accurately recovered, demonstrating that our 
data analysis method is successful at estimating and subtracting 
the stellar contribution. 

Finally, the one-dimensional spectrum of the companion is 
obtained by summing the signal in each spectral channel over 
a four pixel-wide aperture centered on the position of the com- 
panion. The same procedure is applied to the telluric calibrator 
observed after SCR 1845. The companion spectrum is then di- 
vided by the calibrator spectrum and multiplied by a black body 
at 44 500 K (Te/f of the telluric calibrator) to obtain the final 
spectrum that is compared to other methods and used for the 
spectral analysis. 



3.3. Impact of the over-subtraction 

The SD data analysis method relies on two assumptions: (1) the 
data is sufl&ciently well sampled spatially to ensure that it can be 
properly rescaled by interpolation without introducing artifacts 
into the data, and (2) the star "model spectrum" can be properly 
estimated and then fitted in each spatial channel to remove the 
star contribution. If these two conditions are not met, the esti- 
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mation of the stellar halo and speckles will be biased, and the 
companion spectrum will not be recovered properly. 

In our data, the condition (1) is far from being met in J- 
band, where the -54 mas/pixel spatial sampling is 4.2 times 
larger than the Shannon sampling (12.9 mas/pixel). The impact 
of this under- sampling is visible in Fig. 2 (bottom panels) as 
horizontal stripes appearing on both sides of the companion sig- 
nal in /-band owing to the interpolation. With a 27 mas/pixel 
scale, the //-band was undersampled by a factor of only 1.4, 
and no clear indication of under- sampling was visible in the data 
analysis. Although under- sampling is problematic from a signal- 
processing perspective, we verified that it is not the main limi- 
tation here by binning the //-band data to the same spatial sam- 
pling as the /-band data. The resulting spectrum was similar to 
the one without binning within 5%. 

The condition (2) can be met by masking the companion sig- 
nal in the data analysis (as detailed in Vigan et al. 2008). At a 
given angular separation, the size of the mask is limited by the 
value of the bifurcation point (Thatte et al. 2007). The bifurca- 
tion point represents the smallest angular separation at which the 
SD can be applied, which depends on the wavelength range of 
the data and the spatial extent of the companion signal. It mea- 
sures the separation at which the spatial rescaling of the data will 
move the companion signal by at least its spatial extent, allow- 
ing us to fit the stellar speckles and halo on either side of the 
companion signal (see Vigan et al. 2008, Fig. 2 for an example 
of rescaled data). For companions at separations below the bifur- 
cation point, the spatial rescaling will not move the companion 
signal enough to fit the speckles, and the SD will not work. 

If the wavelength range and the separation of the compan- 
ion are fixed, as in our observations, the bifurcation relations 
(Thatte et al. 2007, Eq. 2 and 3) can be used to calculate the 
largest size of the mask that can be applied to the companion 
signal to hide it while still being able to use SD. If the size of the 
mask is larger than the maximum set by the bifurcation point, 
there will not be enough pixels to estimate the star contribution 
and the final result will be biased. Although the mask size was 
chosen to be an adequate size, the eff'ect of the companion sig- 
nal still influences the data analysis outside the masked region. 
A possible explanation is that since NAOS is not a high-order 
AO system, the PSF is not fully diflraction-limited, so it has 
a full width at half-maximum (FWHM) larger than AID. The 
FWHM was estimated to be -100 mas at 1.65 yum in our imag- 
ing data, while the slit has a width of only 86 mas, meaning that 
the companion PSF is larger than the slit. Difl'raction by the edge 
of the slit then broadens the signal on the detector. The impact is 
stronger in /-band because of the poorer seeing, higher airmass, 
and smaller Strehl ratio inducing a larger broadening of the PSF. 

To cope with the over- subtraction, a fake companion spec- 
trum is introduced into the data at a position symmetric to that 
of the real companion with respect to the star. To be as accurate 
as possible, the level of signal of the fake spectrum in each spec- 
tral channel is chosen to match the level of the true companion in 
the same channel (including a stellar contribution). After SD, the 
fake spectrum is extracted using exactly the same procedure as 
described in Sect. 3.2, and the efficiency (i.e. the amount of con- 
served flux) is obtained by dividing the input fake spectrum by 
the extracted spectrum. Finally, the real companion spectrum is 
corrected for the over- subtraction by dividing it by the efficiency. 
The spectrum corrected for the over- subtraction is hereafter re- 
ferred to as the SD-hcompensation spectrum. The spectrum that 
is not corrected will be referred as the SD-only spectrum. A com- 
parison of both spectra in the //-band is plotted in Fig. 5. The 
main changes occur in the spectral slope above 1.6 yum, where 



Table 2. Value of 6rei when comparing the SD-hcompensation 
spectrum to other spectra. 



Comparison spectrum 


/-band 


//-band 


No star subtraction 


28.7% 


17.3% 


Alternative method 


15.7% 


3.9% 


Kasper et al. (2007) 


35.2% 


14.9% 



the variation is limited to a few percent. As we demonstrate in 
the next section, the compensation of the flux loss is necessary 
to match the spectral extraction obtained with other methods. In 
conclusion, we note that our method for compensating for the 
flux loss is valid here because of the modest contrast between 
the companion and the primary, which allows an accurate esti- 
mation of the companion signal. For companions much fainter 
than the star, alternative solutions will need to be developed to 
compensate for the over- subtraction if this occurs. 

3.4. Validation of the method 

To validate the SD, we need to compare the extracted spectrum 
to one obtained after applying an independent method to the 
same data. For that, we use an alternative method based on an 
iterative process: first the stellar halo is estimated and subtracted 
from each spectral channel using a Moffat profile, which closely 
reproduces the shape of the stellar profile; the position of the 
companion is then measured precisely in each channel, and its 
signal is duplicated by symmetry around its center; the dupli- 
cated spectrum of the companion is subtracted from the original 
spectrum to obtain a stellar spectrum that is free of the compan- 
ion contribution; and finally the stellar halo is estimated for this 
new spectrum using a Moffat profile in each channel, and then 
subtracted from the original spectrum. This method is analogous 
to that of Mohanty et al. (2007) presented in Sect. 3.1. As for the 
SD, this method provides a 2D spectrum whose stellar contribu- 
tion has been estimated and removed, but as we see in Sect. 5 
it will be of limited use for very high-contrast in observations 
made by future spectro-imagers. 

The comparison between the different spectra for the /- 
and //-bands is plotted in Figure 6. For reference, the spec- 
trum of the companion for which the stellar contribution was 
not subtracted is also plotted. To quantify the match between the 
SD-h compensation spectrum and other spectra, we used a rela- 
tive spectro-photometric error defined by 



^rel = 



z 



(1) 



where /^^ and Ip are respectively the intensities at wavelength 
p of the SD-h compensation spectrum and the spectrum to which 
we compare it, and N^. the number of wavelengths. Equation 1 
measures the relative difference between the spectral energy dis- 
tributions of the two spectra, up to a scaling constant, making the 
measurement independent of the scaling of both spectra. Table 2 
gives the values of 6rei when comparing the SD-hcompensation 
spectrum to other spectra. We also used a classical coefficient 
to compare the spectra together and found that the trends were 
identical to the 6rei factor, so the values have not been reported 
in Table 2. 

The agreement between the SD-hcompensation and the al- 
ternative method is visually excellent, particularly in //-band 
where the spectral slope and the molecular and atomic features 
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Fig. 6. Comparison of the /-band (left) and //-band (right) spectra of SCR 1845 B obtained without star subtraction (black), with 
SD-Fcompensation (red), and with the alternative data analysis method (blue). The spectrum published by Kasper et al. (2007) is 
also plotted (green). 



are identical between the two methods. This is confirmed by the 
value of 6rei, which is less than 4%. In /-band, the strong over- 
subtraction induced by the SD does not allow us to recover a 
spectrum over the full spectral range. The large over- subtraction 
between 1.10 yum and 1.18 yum causes the signal to drop to neg- 
ative unphysical values. In the properly recovered parts of the 
spectrum, the agreement is good above 1 .2 yum, but below 1 . 1 yum 
the over- subtraction is not completely compensated for. The mis- 
match between the spectra is confirmed by the value of 6rei being 
four times higher than in the //-band, which makes us conclude 
that the /-band SD-hcompensation spectrum is not fully reliable. 

We can then conclude that for data properly sampled and 
taken in good observing conditions, as in the case of our H- 
band data, the SD allows us to attenuate the stellar contri- 
bution and recover the overall spectral shape and features of 
the companion. In particular, we note that the spectrum ex- 
tracted without subtracting the star does not match either the 
SD-hcompensation or the alternative method spectra: the over- 
all spectral slope is diff'erent in both the /- and //-bands, and the 
depth of some strong spectral features in //-band (e.g. the ab- 
sorption at 1.67 yum) is slightly more important after the stellar 
contribution has been removed. This is confirmed in //-band by 
a value for 6rei of 17.3%, which denotes a large discrepancy be- 
tween the SD-hcompensation spectrum and the spectrum without 
stellar subtraction. The use of SD on our data analysis method 
is then necessary to avoid any mis-estimation of the shape of the 
continuum. 



3.5. Comparison to other data 

The spectrum of SCR 1845 B pubHshed by Kasper et al. (2007) 
is also plotted for comparison in Fig. 6. Its diff'erence from the 
SD spectrum is clearly evident, particularly in /-band where the 
H2O absorption between 1.1 yt/m and 1.2 yum is much weaker in 
our data, leading to a value of 35.2% for 6rei. In //-band, the 
visual agreement is much better for both the overall spectral 
slope and the depth of the features, but the value of 6rei is still 
large (14.9%). This discrepancy, mostly due to a small diff'erence 



of slope above 1.6 yum, is compatible with slit losses, a well- 
known issue aff'ecting any AO-assisted long-slit spectrographs 
(Goto et al. 2003). The eff'ect of slit losses generally translates 
into the data as variations in the spectral slope of a few per- 
cent when observing conditions change. Given that Kasper et al. 
(2007) data was taken at a diff'erent date, with a diff'erent instru- 
mental setup, we can then conclude that the agreement with our 
data is reasonable in //-band and does not question the validity 
of the SD method. In /-band, it is difficult to draw any conclu- 
sion as there is a large discrepancy (6rei - 15.7%) between both 
our methods. 

Finally, we note that the data of Kasper et al. (2007) have a 
much lower signal-to-noise ratio (S/N) than our data. They com- 
pleted a 900 s total exposure time on target, but had to use a 
dichroic transmitting only 10% of the incoming light to the sci- 
ence detector to avoid saturation of the primary on the slit, re- 
sulting in an effective exposure time of only 90 s. In contrast, 
our observations totaled 560 s of exposure time with a dichroic 
transmitting 90% of the incoming light to the science camera, re- 
sulting in a much longer effective exposure time of 504 s. While 
this explains the noisier part of their spectrum below 1.2 yum, it 
cannot fully account for the differences between our SD and the 
alternative method spectra, which are probably due to the differ- 
ences in observing conditions. However, it is interesting to point 
out that our observing strategy allows us to use longer exposure 
times at higher throughput, resulting in higher overall S/N. 

4. New spectral analysis of SCR 1845 B 

The significant proper motion of the companion prompted us 
to use the Levenberg-Marquardt algorithm (Press et al. 1992) to 
search for the orbital solution providing the minimal hence 
best-fit to our astrometric measurements. However, we did not 
find any valid orbital solutions. Although our astrometric sam- 
pling remains sparse, additional monitoring will be extremely 
important to confirm that this system is indeed physically bound. 
In the absence of results from astrometry, we then focus the cur- 
rent section on the results obtained from the spectral analysis of 
our new spectrum of SCR 1845 B. 
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Fig. 7. Comparison of the SCR 1845 B /- and if -band spectrum to the synthetic spectra of the BT-SETTLIO library for Teff from 
700 K to 1300 K, log g from 4.5 to 5.0, and [M/H] = dex. The model at Teff = 700 K and log g = 5.5 dex is currently unavailable 
from the library. 
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Table 3. Spectral analysis of our different spectra of SCR 1 845 B 
using the BT-SETTLIO library. 



Table 4. Spectral index values. 



Spectrum 


Method LS^ 


/-band 


//-band 






^eff 


logg 




logg 


SD-hcompensation 


Classic 


900 


4.0 


1100 


4.5 


SD-hcompensation 


Weighted 


1000 


4.0 


900 


5.5 


SD only 


Classic 


600 


3.5 


1000 


5.0 


SD only 


Weighted 


600 


3.5 


800 


5.5 


No star subtraction 


Classic 


1200 


4.0 


1200 


5.0 


No star subtraction 


Weighted 


1200 


4.0 


1200 


5.0 


Alt. method 


Classic 


900 


4.0 


1100 


4.5 


Alt. method 


Weighted 


1100 


4.0 


1100 


5.0 


Kasper et al. (2007) 


Classic 


1000 


4.5 


1100 


5.0 


Kasper et al. (2007) 


Weighted 


900 


4.5 


900 


5.5 



The least squares minimization is either "classic" or "weighted" as 
described in Mohanty et al. (2007). 
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Fig. 8. H2O and CH4 spectral indices defined by Burgasser et al. 
(2006) for the classification of T dwarfs calculated in the spectra 
available for SCR 1845 B. For comparison, the grey dots rep- 
resent the spectral index values for T dwarfs from the SpecX li- 
brary, and the black dots the values for the reference objects from 
Burgasser et al. (2006, Table 4) with available SpecX spectrum. 
The spectral coverage of the Kasper et al. (2007) spectrum does 
not allow us to measure the CH4-/ index. The values of all the 
spectral indices are summarized in Table 4. 



Spectrum 


H2O-/ 


H2O-// 


CH4-/ 


CH4-// 


SD+ compensation 


0.164 


0.252 


0.866 


0.344 


No star subtraction 


0.331 


0.338 


0.557 


0.440 


Alt. method 


0.192 


0.279 


0.480 


0.355 


Kasper et al. (2007) 


0.156 


0.297 




0.283 



ering Teff from 700 K to 1300 K and log g from 4.5 to 5.5 dex, 
and Table 3 (lines labeled "SD+compensation") gives the diff'er- 
ent estimations of Te/f and log g for our SD+compensation spec- 
trum. We also calculated the H2O and CH4 spectral indices de- 
fined by Burgasser et al. (2006) for the classification of T dwarfs 
and present our results in Fig. 8. 

In //-band, the 1.48-1.57 jim and 1.72-1.80 jim regions of 
the companion spectrum are simultaneously reproduced by mod- 
els at Teff = 900-1100 K. The band-head shape between 1.48 
and 1.62 jim is most closely matched at Teff = 900 K. The mod- 
els however fail to reproduce the strength and the details of the 
absorption from 1.62 jim to 1.73 jim. This mismatch is also seen 
in Fig. 23 of King et al. (2010) for 6 Indi Ba/Bb and is inter- 
preted as the result of incomplete CH4 line lists. The mismatch 
is less significant at high surface gravities (log g = 5.0-5.5). As 
a consequence, we present only results when this region of the 
spectrum was ignored in the comparison to the BT-SETTLIO li- 
brary. 

The analysis of the /-band spectrum is more tentative given 
the uncertainties associated with the extraction (Sect. 3.4). The 
models notably confirm a lack of flux between 1.0 jim and 
1.18 fim in the companion spectrum. We note however that the 
shape of the spectrum between 1.0 fim and 1.1 yum can only be 
reproduced by models with Teff > 800 K and log g > 4.5. We 
then finally estimate Teff = 1000 ± 100 K and log g > 4.5 for 
SCR 1845 B. 

The spectral indices calculated in our spectrum (Fig. 8, 
brown line) closely agree with the previously assumed T6 spec- 
tral type. The H2O-/, II2O-H, and CII4-H indices yield a T6 ± 1 
spectral type. The CH4-/ index could not be calculated because 
of the strong over- subtraction in 1 .3 1-1 .34 jim that prevented us 
from recovering a usable spectrum. However, for that spectral 
index, our alternative method, which did not show a flux deple- 
tion as large as the SD spectrum, yielded a T4-T5 spectral type 
that agrees with the other indices. 

The spectroscopic Teff agrees well with the value of 
Teff = 1042 ± 162 K found using the optical/NIR conversion 
scale of Stephens et al. (2009) and assuming a T6 ± 1 spec- 
tral type. Our atmospheric parameters are also consistent with 
those found for field dwarfs by Testi (2009) for this spectral type 
range (1000-1400 K), and based on AMES-COND modelx fit- 
ting (Allard et al. 2001). Finally, they also agree within the error 
bars^ with the values reported by Kasper et al. (2007). 



4.1. Atmospheric parameters 

We compared our extracted SD+compensation spectrum to the 
state-of-the-art BT-SETTLIO library ^ (Aflardetal. 2010) fol- 
lowing the least- square minimization method of Mohanty et al. 
(2007). We restrained the fit to 500 K < Teff < 1500 K and 3.5 
< log g < 5.5, assuming metallicity [M/H] = dex. The best-fit 
solutions were checked visually. Figure 7 presents the compari- 
son of our SD+compensation spectrum to synthetic spectra cov- 

^ http : //phoenix . ens- lyon . f r/Grids/BT- Settl/AGSS2009/; 



4.2. Luminosity, radius, and mass 

We found a luminosity log(L/L0) = -5.28 + 0.08 dex using 
the parallax of Henry et al. (2006), BCk of Golimowski et al. 
(2004) for a T6 + 1 (2.20 + 0.12 mag), the /T-band magnitude 
from Kasper et al. (2007), and their associated errors. The values 
of Teff and luminosity were combined to derive semi-empirical 
radii using the Stefan-Boltzmann formulae L = AnR^crT^rr. With 



^ Note that Kasper et al. (2007) do not state clear error bars for Teff 
OTA)^ g in their work, so our results agree within our error bars. 
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Table 5. Summary of the results obtained with the different spec- 
tra. 





^eff 




SpT 




(K) 


(cgs) 




SD-hcompensation 


900-1100 


4.0-5.5 


T5.0-6.5 


No star subtraction 


1200 


4.0-5.5 


T2.5-4.5 


Alt. method 


900-1100 


4.0-5.0 


T4.0-5.5 


Kasper et al. (2007) 


900-1100 


4.5-5.5 


T5.0-6.0 


= 1000 ± 100 K and 


log(L/Lo) 


= -5.28 


± 0.08 dex, 



estimate a radius R = 0.7 ± 0.1 Rjup. We find a similar value 
(R = 0.7+^-^ Rjup) considering Teff = 1042+162 K. Evolutionary 
models predictions of Baraff'e et al. (2003) are consistent with 
our luminosity, Te/f, and log g estimates for ages > 0.8 Gyr. The 
predicted and the semi-empirical radii agree for ages > 1.5 Gyr. 
These constraints on the system age are fully consistent with the 
one derived by Kasper et al. (2007). Assuming an age between 
1.5 Gyr and 10 Gyr, our values of Te/f and luminosities corre- 
spond to masses between 25 and 65 Mjup. This mass range also 
corresponds to the estimates made by Kasper et al. (2007) using 
Burrows et al. (1997) evolutionary tracks. 

5. Discussion 

Throughout this work, we have demonstrated that the SD 
method is efficient at suppressing the stellar light contribution 
in high-contrast spectra. We now discuss the impact of this sup- 
pression on the determination of the spectral type and the spec- 
tral analysis. We also discuss the advantages of SD compared to 
other methods that have been used in different contexts. 

5.1. Impact of suppressing the stellar contribution 

As we have seen at the end of Sect. 3.4, the SD+compensation 
spectrum does not match the spectrum extracted without sub- 
tracting the star. We quantified it mathematically using the 6rei 
factor (Eq. 1), but it is also important to quantify the impact of 
a biased spectrum on the spectral analysis. For this purpose, we 
performed the spectral analysis (Table 3) and calculated spectral 
indices (Fig. 8) in the different spectra that were available: the 
SD+compensation spectrum, the SD-only spectrum, the spec- 
trum without subtracting the star, the spectrum obtained with 
our alternative method, and finally the spectrum of Kasper et al. 
(2007). The only exception is that spectral indices were not cal- 
culated for the SD-only spectrum because of the unphysical neg- 
ative values that would produce spurious values. 

The Teff (900-1100 K) and log g (4.5-5.5) estimated from 
the SD+compensation and the alternative method are in good 
agreement with the Kasper et al. (2007) results and with the new 
fits we obtained from their data (bottom of Table 3). However, 
their /-band data provides a closer fit to the BT-SETTLIO mod- 
els than our /-band data below 1.18 yum, yielding a value of log g 
that is slightly higher (4.5), which is more in agreement with 
the //-band data. The H2O-/, II2O-H, and CII4-H spectral in- 
dexes are also in excellent agreement (differing by at most 0.05). 
We can then conclude that the overall agreement between our 
SD+compensation spectrum and previously published results is 
good for the estimation of the physical parameters. 

In contrast, the spectral analysis performed on the SD- 
only spectrum and the spectrum without stellar light subtrac- 
tion shows a strong bias. The former yields a very low T^// and 
log g (600 K, 3.5 dex) in /-band and reasonable values in H- 



band, while the latter is clearly biased towards higher Teff in 
both the /- and //-band spectra. For the spectral indices, the 
spectrum with no star subtraction systematically leads to spec- 
tral types lower than T6 (T2-T6). In particular, the H2O-/ index 
is clearly incompatible with the SD+compensation, alternative 
method, and Kasper et al. (2007) spectra. 

Finally, in Table 5 we have reported the estimations of Teff, 
log g, and spectral types (derived from the spectral indices) with 
the different spectra. The spectrum without star subtraction is 
clearly completely biased, while the SD+compensation spec- 
trum is in good agreement with the alternative method spectrum 
and the Kasper et al. (2007) data. 

The conclusion of this comparison is that in high-contrast 
spectroscopic data it is essential to fight two major effects: the 
contamination of the companion spectrum by the star, and the 
systematic flux losses caused either by the slit losses or the a 
posteriori data analysis. Both these effects clearly introduce a 
significant bias in the determination of the spectral type and the 
physical parameters. 

5.2. Comparison of SD to other methods 

In Sect. 3.1, we presented a brief summary of the data analysis 
methods that have been used to differentiate point sources spec- 
tra from a constant or varying background, specifically for high- 
contrast spectroscopy. Given the data they have been applied to, 
these methods are perfectly valid for subtracting the contribu- 
tion of the star at the position of the companion, and they would 
probably provide identical results to our data for SCR 1845 b. 
Our alternative data analysis method is indeed inspired by these 
methods and produces an output similar to previously published 
results. The SD data analysis method is then perfectly compara- 
ble to these other methods in terms of the quality and fidelity of 
the spectral extraction. However, the aforementioned methods 
all rely on the same assumptions: the stellar halo is symmetri- 
cal with respect to the star, and it can be accurately modeled 
using simple mathematical functions (polynomials, Gaussian or 
Moffat profiles, ...). While these assumptions may be reasonable 
for the current state of high-contrast instrumentation, this will 
certainly not be the case for upcoming extremely high-contrast 
spectro-imagers. 

In very high-contrast observations, with extreme AO sys- 
tems and state-of-the-art coronagraphs, the data will be limited 
by the presence of speckles induced by instrumental aberrations. 
This will be the main limitation of any type of data (imaging 
or spectroscopy), in particular because of the temporal variance 
in the speckle pattern (e.g. Hinkley et al. 2007). In LSS data, 
the speckles will modulate the shape of the star halo by creat- 
ing oblique lines of varying intensity along the spectral dimen- 
sion (Sparks & Ford 2002). This is a major limitation because 
the signal of planetary companions will likely be at the level or 
below that of the speckles. In this kind of data, simple data anal- 
ysis methods will not work: although the overall shape of the star 
halo is symmetrical and can be modeled with simple functions, 
the speckles superimposed over the halo are asymmetrical and 
cannot easily be modeled. 

In this context, the SD is the method of choice to remove 
the star contribution (Vigan et al. 2008), be it halo or speckles. 
In very high-contrast observations, using fake companions as in 
Sect. 3.3 to measure the over- subtraction is impossible because 
the signal of the true companion is too faint to be measured ac- 
curately. However, the over- subtraction should be much less sig- 
nificant for several reasons. The main one is that extreme AO 
systems will provide truly diffraction-limited images in the near- 
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IR (Strehl ratio > 80%), hence the companion PSF will be fully 
resolved within the slit. For bright companions, the central peak 
of the companion PSF (which contains -80% of the total en- 
ergy) will then be much easier to mask in the data analysis pro- 
cess and will not introduce any significant bias. The second rea- 
son is that in SPHERE/IRDIS, the slit will be oriented along the 
star-companion axis. This has an impact on the estimation of the 
model spectrum that is fitted in each spatial channel to estimate 
the contribution of the star. In a configuration where the star is 
kept within the slit, the diff'raction pattern and speckles from the 
star remain within the slit when wavelength increases, so they 
have the same contribution at all wavelengths. In the configu- 
ration that was adopted with NACO, these star-related artifacts 
cross the slit because this slit is tilted with respect to the star- 
companion axis, hence their contribution varies much more with 
wavelength: a speckle that is located within the slit at wavelength 
A will move out of the slit at wavelength A -\- AA. This eff'ect is 
small in our data because the variation of the halo is small with 
wavelength, but the presence of bright speckles would bias the 
result even more. 

6. Conclusion 

The spectral characterization of sub- stellar companions is par- 
ticularly important for the understanding of their composition 
and formation processes. The large contrast ratio of these ob- 
jects with their parent star is usually a limitation of our abil- 
ity to extract a clean spectrum, free from any significant con- 
tribution of the star. In the context of the future high-contrast 
specro-imager SPHERE/IRDIS, we have developed a data anal- 
ysis method designed to estimate and remove the contribution of 
the stellar halo and speckles. In the present work, we have been 
able to demonstrate the use of this method with VLT/NACO data 
of SCR 1845 B. 

We have shown that the SD method allows proper recovery 
of the overall spectral slope, as well as the depth of the spec- 
tral features. When good conditions are met (proper sampling, 
good observing conditions), it provides a result identical to an 
independent extraction method and very close to the spectrum 
of SCR 1 845 B previously published. We demonstrated that the 
spectral analysis is strongly biased if the stellar contribution is 
not removed, asserting the essential need for the data analysis 
step. However, the over- subtraction of the companion signal in 
our data makes it necessary to use a fake companion spectrum to 
estimate and compensate for the loss of flux. This efl'ect is partly 
due to the companion PSF being larger than the slit, particularly 
so at lower wavelengths where the AO correction is less efficient. 

From our SD-hcompensation spectrum, we infer a spectral 
type of T6 ± 1, which agrees well with previous studies and 
we estimate Te/f = 1000 ± 100 K and log g = 4.5-5.5 dex by 
comparing our data to the state-of-the-art BT-SETTLIO library 
of models. With a luminosity log(L/Lo) = -5.28 ± 0.08 dex, 
we estimate a radius R = 0.7 Rjup, which is consistent with 
evolutionary models for ages older than 1.5 Gyr 

Finally, these results represent an important milestone for the 
use of the SD method at very high-contrast for the spectral char- 
acterization of planetary-mass objects with LSS and coronag- 
raphy in VLT/IRDIS. Data of future high-contrast imagers and 
spectrographs will indeed be modulated by the presence of in- 
strumental speckles that cannot be estimated and removed using 
simple methods, increasing the need for more sophisticated data 
analysis methods. 
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